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The future finite-time singularities emerging in alternative gravity dark energy models are classi- 
fied and studied in Jordan and Einstein frames. It is shown that such singularity may occur even 
in flat spacetime for the specific choice of the effective potential. The conditions for the avoidance 
of finite-time singularities are presented and discussed. The problem is reduced to the study of a 
' scalar field evolving on an effective potential by using the conformal transformations. Some viable 

' modified gravity models are analyzed in detail and the way to cure singularity is considered by in- 

, troducing the higher-order curvature corrections. These results maybe relevant for the resolution of 

■ the conjectured problem in the relativistic star formation in such modified gravity where finite-time 

' singularity is also manifested. 
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INTRODUCTION 



Several assumptions of the Cosmological Standard Model have been ruled out by the advent of the so called 
precision cosmology capable of probing physics at very large redshifts. The old picture, based upon radiation and 
baryonic matter, has to be revised. Beside the introduction of dark matter, needed to fit the astrophysical dynamics 
■ at galactic and galaxy cluster scales (i.e. to explain clustered structures), a new ingredient is requested in order to 
explain the observed accelerated behavior of the Hubble flow: the so called dark energy. Essentially, data coming 



' from the luminosity distance of la Type Supernovae [1[ , the deep and wide galaxy surveys [.2|| and the anisotropy of 

i Cosmic Microwave Background fl] suggest that the so-called Cosmological Concordance Model (ACDM) is spatially 

fS| ■ flat, dominated by cold dark matter (CDM ~ 25%) and dark energy (A ~ 70%). The flrst ingredient should be able 

' to explain the dynamics of clustered structures while the latter, in the form of an "effective" cosmological constant, 

' should give rise to the late-time accelerated expansion. 

Q\ ' The cosmological constant is the the most relevant candidate to interpret the cosmic acceleration, but, in order to 

1^*) ■ overcome its intrinsic shortcomings associated with the energy scale, several alternative proposals have been suggested 

^ ' (see recent reviews [1,[^,@|): quintessence models, where the cosmic acceleration is generated by means of a scalar fleld, 

• ^ , in a way similar to the early time inflation, acting at large scales and recent epochs frj; models based on exotic fluids 

^ ' like the Chaplygin-gas [S] , or non-perfect fluids ^ ; phantom flelds, based on scalar flelds with anomalous signature in 



H ] the kinetic term [10|, higher dimensional theories (braneworld) All of these models have the common feature to 
introduce new sources into the cosmological dynamics, but, from an "economic" point of view, it would be preferable 
to develop scenarios consistent with observations without invoking extra parameters or components non-testable (up 
to now) at a fundamental level. 

Alternative theories of gravity, which extend in some way General Relativity (GR), allows to pursue this different 
approach (no further unknown sources) giving rise to suitable cosmological models where a late-time accelerated 
expansion is naturally realized. 

The idea that Einstein gravity should be extended or corrected at large scales (infrared limit) or at high energies 
(ultraviolet limit) is suggested by several theoretical and observational aspects. Quantum fleld theories in curved 
spacetime, as well as the low-energy limit of string theory, both imply semi-classical effective Lagrangians containing 
higher-order curvature invariants or scalar-tensor terms. In addition, GR has been deflnitely tested only at Solar 
System scales while it may show several shortcomings, if checked at higher energies or larger scales. Besides, in the 
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opinion of several authors, the Solar System experiments are not so conclusive to state that the only reliable theory 
at these scales is GR. 

Of course modifying the gravitational action asks for several fundamental challenges. These models can exhibit 
instabilities or ghost -like behaviors while, on the other side, they should be matched with observations 
and experiments in the low energy limit (in other words. Solar System tests and Parameterized Post Newtonian 
(PPN) limit should reproduce the results of GR appropriately). Despite of all these issues, in the last years, several 
interesting results have been achieved in the framework of the so called modified gravity at cosmological, galactic and 
solar system scales (see Refs. EBl for review). 

For example, there exist cosmological solutions that lead to the accelerated expansion of the Universe at late times 
in specific models of f(R) gravity as is discovered in refs.[ll, [l3, In some of realistic theories of this sort the 
problems indicated in [3| maybe overcomedjl^. 

There exist viable f{R) models that can satisfy both background cosmological constraints and stability conditions 
[20l . [2ll . [2^ . [23 . [23 . IsqI . [iol as well as local tests. Recently many works have been devoted to place constraints 
on /(^)-models using the observations of Cosmic Microwave Background anisotropies and galaxy power spectrum 

Besides, considering /(i?)-gravity in the low energy limit, it is possible to obtain corrected gravitational potentials 
capable of explaining the flat rotation curves of spiral galaxies and galaxy cluster haloes without considering huge 
amounts of dark matter jlS, HI, [13, S, lU, [H, |3§| and, furthermore, this seems the only self-consistent way to 
reproduce the universal rotation curve of spiral galaxies |33|. On the other hand, several anomalies in Solar System 
experiments could be framed and addressed in this picture |34|- 

However, a fundamental task which has to be faced for any alternative gravity model is to classify singularities 
which could emerge at finite time and propose the way to avoid it. From a physical point of view, this point is crucial 
in order to achieve viable and self-consistent models, especially in the possible applications. 

In this paper, we discuss the future singularities which can, in principle, appear in dark energy models coming from 
alternative gravity theories (higher-order or non- minimally coupled gravity). 

In fact, when dark energy models with the effective equation of state parameter close to —1 were added to the 
list of admissible cosmological theories to explain the observed accelerated behavior, due to to violation of all (or 
part) of energy conditions, strange features emerged in the future. For instance, it is well-known that phantom dark 
energy brings the universe to finite-time Big-Rip singularity fss','^?!. Moreover, the effective quintessence dark energy 
cosmologies may end up in (softer) finite-time singularity [38, 41]. For such effective quintessence dark energy models 
only part of energy conditions does not work in the standard way. Nevertheless, they show up the rip singularity 
behaviors which have been classified in Ref.fs^l- 

It is clear that, qualitatively, the same situation should occur also in modified gravity cosmologies Indeed, it is 
quite well-known that some versions of modified gravity (Hke f{R)) have an effective ideal fiuid description jH]. Hence, 
precisely the same singular behavior should be typical for the (effective phantom/quintessence) modified gravities in 
future too. Indeed, it was found some time ago [3g| that modified gravity becomes invalid (complex theory) at the 
point where mathematically-equivalent scalar-tensor dark energy theory enters to the Big Rip singularity. Moreover, 
the effective phantom behavior may enter a transient phase and future singularity does not occur if some higher order 
terms (like R^) are added to initially phantom- like models [i^. The same approach has been recently considered in 
ref.fi^l to remove the singularity in order to avoid the conjectured problems with neutron stars formation in modified 
gravity. 

In this paper, we want to discuss, in general, the problem of finite-time singularities and discuss some ways to avoid 
them in viable models which well-fit data at local and cosmological scales. 

The layout of the paper is the following. In SectUll we describe, in general, the problem of finite-time singularities 
in dark energy models coming from f{R) to scalar-tensor modified gravities. Sect. Hill is devoted to the conditions 
for singularity avoidance in /(i?)-gravity and in its scalar-tensor counterpart. In Sect. IIV^ we discuss the singularity 
problem in some physically viable /(i?)-models adopting a conformal transformation approach. This method allows, 
in principle, to discriminate singularities by studying the behavior of effective scalar field potential after dynamics 
has been conformally reduced to the Einstein frame. In particular, we study the effect of adding a correction term, 
proportional to i?" with rt > 2, to modify the structure of the potential at large values of R and cure the singularity. 
Discussion and conclusions are drawn in Sect. [Vl 
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II. FINITE-TIME SINGULARITIES IN DARK ENERGY MODELS: FROM F(i?)-GRAVITY TO 

SCALAR-TENSOR THEORY 



Let us start with a generic action of F(i?)-gravity which is a straightforward extension of General Relativity: 



d'^V-9[^+C^\ ■ (1) 



Here F{R) = R + f{R) is an appropriate function of the scalar curvature R and £m is the Lagrangian density of 
matter. By the variation over the metric tensor g^,y, we obtain the fourth-order field equations: 

1 

-9^uF{R) - R^,F'{R) - g^,,UF'{R) + W^W,F'{R) = -yT(„)^, . (2) 

In Eqs.([2]), T'(m);ijy is the matter energy-momentum tensor. Contracting Eqs. ([2]) with respect to /j, and v, we obtain 
the trace equation: 

^2 

2F(R) ^ RF'{R) -iUF'iR) ^ -—T . (3) 
To recover, formally, General Relativity, the Eq.|[3l) can be rewritten as 

i? + 2/(i?)-i?/'(i?)-3n/'(i?) = -yr. (4) 

In order to study how finite-time singularities emerge and can be discussed, let us consider, for the moment, classes 
of models which are paradigmatic for our purposes. For example, in ref.|24|, a model which easily passes local tests 
and several cosmological bounds, has been proposed: 

or otherwise written as 

(—) 

fHsiR) - -XRcy^ (6) 

(^) +^ 

Here m is a proper scale and ci, C2, n and A are dimensionless positive constants {n is not restricted to be an 
integer) and Rc is positive constant. When the curvature is sufficiently large at dark energy epoch, this model can be 
approximated as follows 

fiR)^-2A+^. (7) 



In case of one may identify 



Then Eq.Q reduces to 



m^Ci m^™+^ci , , 

2A = , a = 5 — - . (8) 

C2 ci 



i? + 3an(i?-"-i) -0 . (9) 



Since the large curvature regime is considered, the cosmological constant term appears as a next-to-leading order 
correction, compared with the first term in ([9]), and we neglect it. If we define x 

x=i?-"-^ (10) 

and the FRW metric with flat spatial part is chosen, 

ds'^ ^ -^dt'^ + a{tf {dx'f . (11) 

i=l,2,3 
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Eq.ljQ]) has the following form: 



da 



(12) 



Note X = corresponds to the curvature singularity R oo. Note that as other dark energy models with EoS 
parameter around —1, the above gravitational alternative for dark energy also has the singularity as it will be explained 
below. The fact that such F(i?)-gravity may show the phantom- like behavior has been established in ref.fi^. It was 
demonstrated there that account of (or similar nature term) makes the phantom phase to be transient and removes 
the singularity. In principle, the phantom phase in F(i?)-gravity may end up as Big Rip-like type singularity fH] as 
it was demonstrated in ref.|36|- 

First we consider the classical equation of motion: 



3a 



(13) 

The difference between the cosmological equation lfT2|) and the classical equation ifTSjl is the second term depending 
on H. This second term gives the only sub-leading contribution, which will be shown in the analysis from (flSl) . where 
the H dependence will be explicitly included and it will be shown that the result from the classical analysis here will 
be reproduced. For Eq.lfTS]). one gets an exact solution: 



(14) 



Here C and to are constants. Note 2 > 2(n+l)/(n + 2) > 1. Then x vanishes in a finite time t = to, which corresponds 
to the curvature singularity in (fT3l) . 

We now investigate the asymptotic solution when the curvature is large, that is, x is small. As there is a curvature 
singularity, one may assume 



H^ho (to-ty 



(15) 



Here Hq and /3 are constants, Hq is assumed to be positive and t < io as it should be for the expanding universe. 
Even for non-integer (3 < 0, some derivative of H and therefore the curvature becomes singular. The case /3 = 1 
corresponds to the Big Rip singularity Fig. [TJ Furthermore since (3 — corresponds to the de Sitter space, which has 
no singularity, it is assumed /9 0. When /? > 1, the scalar curvature R behaves as 



Figure 1: Plot of H ^ ho (to — t) ^ . The Big Rip singularity occurs for /3 = 1 and t — to = &. 



R - 12H^ - I2hl (to - t) ^'^ 
On the other hand, when /3 < 1, the scalar curvature R behaves as 

R^6H ^ -6/io/3 (to - ty'^-^ 



(16) 



(17) 



Then R diverges when /3 > -1 but /? ^ 0. 

We now consider three cases (compare with [s^l): 1) /? = 1, 2) /? > 1, 3) 1 > /9 > 0, and 4) > /3 > -1 see Fig. [2l 
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R(t) in case 1 

- R(t) in case 2 
R(t} in case 3 

- R(t) in case 4 



Figure 2: Behavior of R{t) for the four cases in the text. The dash-dot line represents 7? for /3 = 1 in Eq. lfTSll : the dashed line 
is for /3 > 1 in Eq.^; the dot line is R{t) for 1 > /3 > 0; finally, the solid line R{t) is for > /3 > -1 in Eq.^T^ respectively. 



• In case 1) /3 = 1, since 



and therefore, from lfTO|). we find 



and the l.h.s. of (fT2l) behaves as 



„ Uhl + 6ho 

It ^ 7i . 

X + 3Hx^ {to ~ t f 



(18) 



(19) 



(20) 



but the r.h.s. behaves as 



6a 

virhich is inconsistent since the powers of the both sides do not coincides with each other. Therefore, /? ^ 1. 
• In case 2) /3 > 1, we find 

R = 12H^ + 6H ^ 12H^ ~ (to - ty^^ , 

and therefore 

In the l.h.s. of lfT2|) . the second term dominates and the l.h.s. behaves as 

X + 3Hx-3Hx-{to-tf^'"^'^-' . 
On the other hand, the r.h.s. behaves as 

-2/3 



(21) 



^X-^/("+i)^(io-t)- . 
3a 

Then by comparing the powers of the both sides, one gets 

/3(2n + 1) - 1 = -2/3 , 
which gives /3 = l/(2n + 3) but this confiicts with the assumption /3 > 1. 



(22) 
(23) 
(24) 

(25) 
(26) 
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• In case 3) 1 > /? > or case 4) > /3 > -1, we find 

R = 12H^ + 6H r^6H (<o - ty'^~'^ , (27) 

and therefore 

. (28) 
Then in the l.h.s. of ifT2|) . the first term dominates and the l.h.s. behaves as 

^ + 3Hx-X-{to-tf^"^'^+''-' ■ (29) 
On the other hand, the r.h.s. behaves as 

^X-'^^'-^'^-ito-ty'^' ■ (30) 
Then by comparing the powers of the left-hand side and the right-hand side, the consistency gives 

/3(n + 1) + n - 1 = -/3 - 1 or /? = -n/{n + 2) . (31) 

This confiicts with the case 3) < /3 < 1 but is consistent with the case 4) > /? > — 1. In fact, by substituting 
((311) into ([28]), we get 



X^ito- t)2("+i)/(«+2) . (32) 

which corresponds to (fT4| . Since > /? > — 1, this singularity corresponds to Type II in [3^ . 

Thus Jthe sudden finite-time curvature singularity really appears in the Hu-Sawicki model. 

In [38|, there was suggested the classification of the finite-time singularities in the following way: 

• Type I ("Big Rip") : For i ^ t^, a — > od, p — > oo and \p\ — > oo. This also includes the case of p, p being finite at 

• Type II ("sudden") : For t ^ ts, a as, p ^ ps and \p\ oo 

• Type III : For t ^ ts, a as, p ^ oo and \p\ — > oo 

• Type IV : For t — > tg, a Os, p ^ 0, |p| and higher derivatives of H diverge. This also includes the case 
when p (p) or both of them tend to some finite values while higher derivatives of H diverge. 

Here ts, as and ps are constants with Os ?^ 0. We now identify ts with to. The Type I corresponds to /? > 1 or /? = 1 
case and we have a Big Rip singularity [47|, whereas the Type II to — 1 < /3 < case and corresponds to the sudden 
future singularity. Type III to < /? < 1 case and is different from the sudden future singularity in the sense that p 
diverges, and Type IV to /3 < — 1 but (3 is not any integer case. 

Let us now remind that the Type II singularity has been already discussed in several dark energy models [4l| besides 
i^(i?)-gravity. Here we consider several theories in the FRW space-time with fiat spatial part (fTTI) . For the Type II 
singularity the Hubble rate H = a/a has the following form 

H = Ho + Hiitn-ty . (33) 

Here Hq, Hi, to and 7 are constants. We now choose < 7 < 1. Then H is finite H Hq in the limit oi t ^ to but 
H diverges as 

H^m^{to-ty-' , (34) 

which generates the singularity in the scalar curvature R since (see Fig.|[3l)) 

R = 12H^ + 6ij - 6iJi7 {to - ty~'^ . (35) 
We should note that the energy density p is finite since the first FRW equation gives 
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P = , (36) 



Figure 3: Behavior of R as given in Eq. I|35p . We have assumed 7 = |- 



and therefore p Hq < oo in the limit t Iq. Hence, the curvature singularity could occur even if the energy 

density is finite as in some other quintessence models. 

We now give an explicit example of the ideal fiuid which gives the singularity in ([SS)) . First we should note that 
the second FRW equation has the following form: 



p= -\ [2H + 3H 



For the Hubble rate in ||33]), Eqs.JMl and ^ give 

p^^{Ho + H^ {to - tyf 

Then by deleting t in the two equations of ((38| . we find 



p= \Ho + Hi 



P + P = 



{P + P) 



(to - 1) 



7-1 



7/(7-l)N 



(37) 



(38) 



(39) 



Eq. l|39p can be regarded as an equation of state (EoS). Conversely if we consider the perfect fiuid satisfying the EoS 
(|39l) . the singularity (|33|) occurs. 

We now consider the occurrence of singularity ([331 , in terms of the scalar-tensor theory, whose action is given by 



(40) 



where V{(f>) is the scalar potential and is the kinetic function, respectively. Note that for convenience the kinetic 
factor is introduced. We should note that scalar field may be always redefined so that kinetic function is absorbed. 
In the spatially fiat FRW space-time ifTTj) , the energy density and the pressure of the scalar field is given by 



1 



1 



Combining the FRW equations ([SG)) and ([37l) with (|4T|) . one obtains 



y(0) 



1 



We now consider the theory in which V{(t)) and a;(0) are given by 



3H^ + H 



a;(0) 



r/'(0) 



F(0) 



1 



3/(0)' + /'(0) 



(41) 



(42) 



(43) 
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where f{(f>) is a proper function of (j). Then the following solution is found 

<j>^t, Hit) - fit) . (44) 

In case of (l33l). we find 



c^(<^) = -^(io-0r~\ Vi<j>) = ^ i3Ho + 3H,ito-^)^)' + Haito~^)^~' ■ (45) 



Thus, with such potential choice, the singularity l(33|) could be easily realized and the energy density is, of course, finite. 
It is easy to construct the models showing the finite future singularity in other models, say, the scalar-Gauss-Bonnet 
theory, non-minimal theoriesfi^, etc. 

We now consider what kind of scalar potential can generate the singularity in the scalar field in the fiat space-time 
background. As a form of the singularity, we now assume 

^ ^ (f>o + (f>i ito - t)^ . (46) 

Here (j>i, and 7 are constants. One may take 7 to be positive but not integer. Then some derivative of the scalar 
field has singularity. 

In the next section, we will rewrite the i^(i?)-gravity in the scalar-tensor form, where the metric is rescaled as given 
in ((53l) . Due to the scale transformation, the curvature singularity in the original frame of the F(i?)-gravity does not 
appear in the rescaled metric in the the scalar tensor frame. The singularity occurs in the scalar field as follows from 
(|46|) . Hence, such scalar field singularity l(46|) is called finite-time singularity even in fiat Minkowski space. 

In the fiat space-time background, the equation of the scalar field is given by 

4> + V'i(P) = 0. (47) 
Here 1^(0) is a potential of the scalar field. By substituting (|46l) into l(47|) . we find 

(/.i7(7-l)(io-i)^"' + V^'(0) =0 . (48) 

Since l(46|) can be rewritten as 

Do 



to-t={ ] , (49) 

by substituting the expression (|49|) into (|48l) . we find 



</'i7(7-l)(^V^)' '^Vr(0) = O, (50) 



which gives the form of the potential as 



Vi^) = Vo-^(^-^y~'\ (51) 



Hence, it is found that finite-time singularity l(46|) can be reahzed by the potential ([511) even in fiat space-time. Thus, 
we demonstrated that for variety of dark energy models including modified gravity the finite-time singularity easily 
occurs even in the situation when effective EoS parameter is bigger than —1 (the effective quintessence). 



III. THE AVOIDANCE OF FINITE-TIME SINGULARITY IN MODIFIED J^(i?)-GRAVITY 

Let us consider the action ^ again. By introducing the auxihary field A, we rewrite the action |(T|) of the FiR)- 
gravity in the following form: 

[d^xV^{F'iA)iR-A)+FiA)} . (52) 



2k 
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Here we neglect the contribution from the matter. By the variation over A, one obtains A = R. Substituting A = R 
into the action l(52|) . one can reproduce the action in Jl]). Furthermore, we rescale the metric in the following way 
(conformal transformation) [l^ : 

g^.^e^g^., a^~\nF'{A). (53) 
Hence, the Einstein frame action is obtained: 



Via) ^ e^g (e"^) - e^/ {g (e"^)) = ^ ^ |^ (54) 

Here g (e~'^) is given by solving the equation cr = — In (1 + f'{A)) = \nF'{A) as A = g (e""^). In terms of f{R) the 
potential V{(t) could be rewritten as 

(l + f'(A)f 

For the class of models with f{R) behaving as in |[7|), we find 

n-)--2A-^, (56) 



for large scalar curvature. When the curvature is large in the model ([7]), one gets 

an 



By combining l(56|) and (|57| . the potential V{a) is given in terms of a as: 



(57) 



y(cr) - -2A- (n + l)a — . (58) 

Van/ 



By comparing (|56l) with ijSTj) and identifying a with (0 — 0o)/0o, we find 



= 2 - - , (59) 



n + 1 7 
or 

2(n+l) , , 

Since 7 is fractional in general, the scalar field a generates the singularity in l|46p . which now corresponds to the 
curvature singularity in the Jordan frame. We should note that the curvature in the Einstein frame is in general not 
singular. 

Let us consider the realistic models which unify the early-time infiation and late-time acceleration and which were 
introduced in refs. [sol. !40|. 

In order to construct such models, the following conditions are used: 



Condition that infiation occurs: 



lim f{R) = -A, . (61) 

R — ^00 



Here is an effective early-time cosmological constant. 
• The condition that there is fiat space-time solution is given as 

/(O) = (62) 
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• The condition that late-time acceleration occurs should be 

/(i?o) = -2Ro , f'iRo) - . (63) 

Here Rq is the current curvature of the universe and we assume Rq > Rq. Due to the condition (|63l) . f{R) 
becomes almost constant in the present universe and plays the role of the effective small cosmological constant: 
Aeff ~ -/(i?o) = 2i?o. 



An example which satisfies the conditions (|6T|) . l(62|) . and ((63|) is given by the following action[3 

(i?-flo)'"+V<'+^ 
/o + /i{(i?-i?o)'"+V<+^} 



/ D 7? \2ri+l n2ri+l 



Here n is a positive integer. 

The conditions l(6T|) and l|63p require 



f I f r,2n+l " ' f 

Jo + JlRo h 



2i?o , ^ = . (65) 



One can show that for the above class of models, there does not occur the singularity or it is difficult that the 
singularity occurs in visible future. We now work in the scalar-tensor form or the Einstein frame l(54|) . 

In case of the model in i:2l|, in the limit oi A = R ^ -l-oo, f{A) becomes a small constant — Aoff corresponding to 
the small effective cosmological constant in the present universe. In the limit, we find f'{A) — + and the potential 
(|55|) behaves as in fM)) and therefore 

V{a) ^ Aeff . (66) 

Then the value of V is very small and the curvature singularity could be easily generated. In the models exhibiting 
this kind of singularity due to fact that singularity is at finite distance with current energy scale, there was conjectured 
[i^ l that neutron stars cannot be formed. (This is really conjecture because the derivation of stars formation follows 
the same approximation as in the usual Einstein gravity, actually neglecting the higher derivatives terms typical in 
i^(i?)-gravity. This derivation should be reconsidered within the real F(_R)-gravity with account of higher derivative 
terms as well as non-lineari^what is technically not easy task. For instance, the non-linear structure of modified 
gravity leads to oscillations [Hi already in the very simple approximation.) The conjecture is that if the density of 
the matter becomes finite but large and reaches a critical value, the curvature singularity occurs and the star could 
collapse. Therefore the star with density larger than the critical value could not be formed. As is mentioned above 
these simple considerations may be not valid with account of higher derivative terms where star formation process 
should be reconsidered from the very beginning. 

On the other hand, for the class of models satisfying the conditions l|61H63p like the model (|65l) . the singularity does 
not occur easily. For this class of the models, f'{R) vanishes aX R = Rq and in the limit of i? ^ oo. When R = Rq, 
the condition IpT]) gives 

V{a) 2Ro , (67) 

which could be quite small. Since the value of R is finite, l(67|) does not correspond to any singularity. On the other 
hand, in the limit of i? — > -l-oo, by using the condition (|63| . we find 

Via) -> A, . (68) 

Since Ai corresponds to the effective cosmological constant during the inflation, the energy scale is not small, typically 
it is the Grand Unification scale. Therefore the value of V{a) could be very large. Then the singularity could be 
generated only at the energy density larger then the energy density corresponding to the infiation of the early universe 
but it does not occur around the energy-density which is typical for neutron star. 

Even for the class of dark energy models where singularity occurs at smaller energies, there is scenario to avoid the 
singularity proposed in ref.ji^. Indeed, let us consider the model where ,f{R) behaves as 

f{R) ^ /oi?" , (69) 
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with constants /o and a > 1. If the ideal fluid, which could be the matter with the constant EoS parameter w: 
p = wp, couples with the gravity, when the /(i?)-term dominates compared with the Einstein-Hilbert term, an exact 
solution is [4^ 

h , 2a 

a — aot ° , ho = 



3(1 

6/n/in / , 



' 6ho + I2hl)" ' {(1 - 2a) (1 - a) - (2 - a)ho} 



Po 

2 

3(l+«;) 



(70) 



When a — 1, the result ho — 3(1^^1 in the Einstein gravity is reproduced. The effective EoS parameter Wcs may be 
defined by 

3(1 + Wcs) 



By using l|70p . one finds 



l + w , , 

Wcff = -1 + . (72) 



Hence, if w is greater than —1 (effective quintessence or even usual ideal fiuid with positive w), when a is negative, 
we obtain the effective phantom phase where Wcs is less than —1. This is different from the case of pure modified 
gravity. On the other hand, when a > w + 1 (it can be even positive), WcS could be negative for negative w. Hence, 
it follows that modified gravity minimally coupled with usual (or quintessence) matter may reproduce quintessence 
(or phantom) evolution phase for dark energy universe in an easier way than without such coupling. 

If we choose a to be negative in ijGQ]). when the curvature is small, the f{R) term becomes dominant compared 
with the Einstein-Hilbert term. Then from l(72|) . we have an effective phantom even if w > —1. Usually the phantom 
generates the Big Rip singularity. However, near the Big Rip singularity, the curvature becomes large and the Einstein- 
Hilbert term becomes dominant. In this case, we have WcS ~ w > — 1, which prevents the Big Rip singularity. 

One may add extra term to f{R) ([Ml) as (45l.l4i|. 

f{R) = foR" + fiR^ , (73) 

Here we choose /i > and /? > 1. Then for the large curvature, the second term could dominate. Then for the large 
curvature, the potential (|55l) behaves as 

(74) 

Then if 1 < /3 < 2, the potential is positive and diverges near the curvature singularity R ^ 00, which could 
prevent the curvature singularity even if the singularity is the Big Rip type (Type I) or Type Il[^ or other softer 
singularity 01 . 

Thus, we demonstrated that for large class of viable F(i?)-gravities the finite-time singularity occurs in so distant 
future that it cannot infiuence the current universe processes. From another side, there exists the trick introduced in 
refs.^45. i46] how by adding extra term to modified gravity to remove the singularity. 



A. Singularity avoidance in integral modified gravity models 

As it is clear from l(52|) . for F'{R) = 1 + f'{R) > 0, the square of the effective gravitational coupHng is positive. 
However, due to n^g = k^/F'{A), the theory could enter an anti-gravity regime. In order to avoid this anti-gravity 
problem from the very beginning, we may consider a model given by 

fiR) = -/o / dRc ^ . (75) 

Jo 

Here a, (3, fo, and i?i are constants and we assume < i?i ^ A^. In this model, the correction to the Newton law is 
very small. Then by construction, as long as < /o < 1, we find f'{R) > — 1 or F'{R) > and therefore there is no 
anti-gravity problem. Since 

f{Ri) --fa f ' di?e"(«-«i)^" = -/oA„(a)i?i , A„(a) = / dxe"^ , (76) 
Jo Jo 
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and —f{Ri) could be identified with the effective cosmological constant 2Ro, we find 

/oA„(a)i?i = Ro . (77) 
Note that An{0) = 1, An{+oo) — 0, and A'{x) < 0. On the other hand, since 

/(+oo) ~ / dRc'^ = -/o/3A, , (78) 
Jo 

and — /(+oo) could be identified with the effective cosmological constant at the infiationary epoch, A^, we find 

foP=l- (79) 

Then the conditions l(6T|) and (|63| are satisfied. The condition l(6T|) is, of course, satisfied by construction (|75| . As 
discussed around Eq. l|68p . as long as A^ could be large enough, the potential becomes large when the scalar curvature 
R is large and there could not occur or could be difficult to realize the singularity. We also note that if we add the 
second term fiR^ with 1 < /? < 2 in iffS]). the singularity is completely removed since the potential diverges in the 
limit of i? ^ oo. Thus, for large class of viable modified gravity models the singularity maybe easily removed by 
adding of extra term which is actually relevant at the early universe. Moreover, the typical energy scales of neutron 
star and singularity formation process (above the infiation scale) are at large distances and they are not relevant for 
each other. 



IV. THE SINGULARITY PROBLEM IN /(i?)-VIABLE-MODELS 

In this section, we are going to discuss the curvature singularity problem that affects several infrared-modified f{R) 
modelji^]. As we have said above, considering f{R) ^ i?-gravity means that a new scalar degree of freedom has to 
be taken into account. Conformal transformations of the metric can be used to make it explicit in the action [4^.[49|. 

For our goals, we consider a class of /(i?)-models which do not contain cosmological constant and are expHcitly 
designed to satisfy cosmological and Solar-System constraints. In practice, we choose a class of functional forms 
of f{R) capable of matching, in principle, observational data [H3|- First of all, any viable cosmological model have 
to reproduce the CMBR constraints in the high-redshift regime. Secondly, it should give rise to an accelerated 
expansion, at low redshift, according to the ACDM model. Thirdly, these models should give rise to a large mass for 
the scalaron jH] in the high-density region where local gravity experiments are carried out. In such a regime, the 
perturbation in R can be larger than the background value i?o, which means that the linear approximation to derive 
the Newtonian effective gravitational constant in a spherically symmetric spacetime ceases to be valid [H^ . IH^ I . In the 
non-linear regime with a heavy scalaron mass, however, it is known that a spherically symmetric body has a thin-shell 
[13, m, [13, EE] through so called the chameleon mechanism jH, (see also Ref. jH, [6l|). When a thin-shell is 
formed, an effective coupling that mediates the fifth force gets smaller. This allows the possibility that the f{R) 
models which have a large scalaron mass in the high-density region can be compatible with local gravity experiments. 
Then there should be sufficient degrees of freedom in the parameterization to encompass low redshift phenomena (e.g. 
the large scale structure) according to the observations [51|. Finally, small deviations from GR should be consistent 
with Solar System tests. All these requirements suggest that we have to assume the limits 

lim f{R) ~ constant, (80) 

R — >oo 

^im/(i?)-0, (81) 
which are satisfied by a general class of broken power law models, as those proposed in jl^l 



fHs{R) — — Ai?c 



2n 



(82) 



Since f{R = 0) = 0, the cosmological constant has to disappear in a fiat spacetime. The parameters {n, A, Rc} are 
constants which should be determined by experimental bounds. 
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Other interesting models with similar features have been studied in |23l.l24l.l39l.l40l|. In all these models, a de-Sitter 
stability point, responsible for the late-time acceleration, exists for R = Ri (> 0), where i?i is derived by solving the 
equation Rif,n{Ri) = 2/(i?i) [Hi|. 



In the region R^ Rc, model ([821) behaves as 



fhybridl{R) — ~ ^Rc ^ — {Rc/ R) 



(83) 



where s is a positive constant. The model approaches ACDM in the limit R/Rc 
Finally, let also consider the class of models [13, lU, HB] 



fhybrtd2[R) — — ^Rc j ' ^^^^ 

Also in this case A, q and Rc are positive constants (note that n, s and q have to converge toward the same values to 
match the observations). We do not consider the models whit negative q, because they suffer for instability problems 
associated with negative /,_r_r [25| . 



A. The role of conformal transformations 



Conformal transformations can play a key role to classify singularities in modified gravity (apart from general 
classification presented in second section). As we have shown in above section, at scale Rc and beyond, the expansion 
rate of the Universe is set primarily by the matter density, just as in standard cosmology, with small corrections. 
Once the local curvature drops below Rc, according to the chameleon mechanism, the expansion rate feels the effect 
of modified gravity. The spacetime curvature, on the other hand, is controlled by the further scalar degree of freedom 
a which gravity acquires. It obeys the usual scalar field equation with potential V{a), the shape of which is directly 
determined by function f{R), and a driving term from the trace of matter stress-energy tensor. But a problem comes 
out at this point: it turns out that precisely those functions f[R) that lead to Einstein-Hke gravity action in the 
large curvature regime, yield a potential V with an unprotected curvature singularity. (Note that just just the same 
occurs for number of realistic quintessence dark enrgies.) Let us consider, for example, the model (|82l) which has been 
constructed to avoid linear instabilities. The conformal transformation gives a = — ln(l + f'(R)) = — hiF'{R) with 
the potential defined in Eq. ljSS]) . For the model (|82)) . the scalar field a is given by 

a{x) = -\n 2- ^'"'^'"''M (85) 

where the cross-over curvature scale Rc can be reabsorbed into a rescahng of coordinate (which can be assumed 
dimensionless and can be measured in length units corresponding to Rc). Suitable coordinates are x = ^ > and 
Rc ^ Pg ^ 10^^'^g/cm'^ for the Galactic density in the Solar vicinity and Rc ^ Pg ^ lO^^^g/cm'^ for the present 
cosmological density. For the large curvature Hmit R ±oo, we have that a = —0.70 while, for i? ^ 0, it corresponds 
to fj = —0.70, which gives us a hint that the potential is going to be a multi- valued function. For R ^ 1, we have 
a — 0. We have V{a) — > oo for i? — > 1, then it is singular for this curvature value, while for i? ^ oo, we have V{a) — 2 
and fiat spacetime in the limit R ^ 0. The analytic expression is 



(2(a;-nA)a;2" +a;)^ 

For our goals, it is important to study the trend of the potential and of the field a through a parametric plot of two 
functions with respect to x — In this way, singularities can be easily identified (see Fig. 
For the hybrid-1 model, a assumes the form 



^ \ 2s+l 



cr = -In 2-2s - A . (87) 



>From FiglU for A = 2 and s = 1, we can see that we have tr = in the limit R — oo and cr ^ oo for R — 0. 
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Figure 5: Plot of Eq.lEZl) for A = 2 and s = 1. 



The effective potential is 

(^2a _ 2s _ 1) A 



V{a) = - 



(a;2«+i - 2s\f 



and, from a rapid inspection of Fig. [6] with A = 2 and s = 1, we have V — 2 for R ±oo. 
Finally, let us analyze the model l(84|). The scalar field assumes the form 



In + 2 ) (89) 



and then a — oo for R ^ 0. 
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- hybridl potential 



1 1.5 



Figure 6: Plot of Eq.((88l) for A = 2 and s = 1. 



- tiybrid2 scalar field 





x)R/R 



Figure 7: Plot of Eq.JHS]) for A = 2 and <j = 1. 



The potential is 



Via) 



{q + l)xi+^X 



(90) 



which gives V = for R ±oo and for R — 0. As concluding remark, we can say that the conformal transformations 
allow to recast the singularity problem in terms of the scalar field and its potential. This fact suggests an easier 
interpretation of the singularities, at least, in some specific cases. 
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Figure 8: Plot of Ea. l|90p for A = 2 and g = 1, it is worth of noting the double sigmoid trend. 
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Figure 9: Parametric plot of Eqs.(l90|) and for A = 2 and q = 1. 



B. Singularities vs Chameleon mechanism 



In order to compare these results, we can take into account another approach where the matter coupHng is considered 
as also been studied already from Sami et al. [gOl- This allows to compare singularities with density scale and could 
be particularly useful to construct models where finite singularities are avoided at infrared scales. 

Starting again with action (U) which leads to the equations of motion the evolution of the scalar degree of 
freedom is given by the trace: 



2F{R) - RF'{R) - 3DF'{R) = 

Another convenient way to define the scalar function a is 

cr = F\R) ~ 1, 



(91) 



(92) 
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which is expressed through Ricci scalar once F{R) is specified. We can write the trace equation (equation ([3])) in the 
term of V and T as 

^ dV 



The potential can be evaluated using the following relation 

then V{a) is given by the pair of functions {<j{R), V{R)}. Let us consider now the model l(82|) as also is showed in 
[60| . The scalar field a is given by 

We can compute V{R) for a given value of n. In the case of n = 1 and A = 1, we have 

where x = R/Rc- In the FRW background, the trace equation can be rewritten in the convenient form 

dV k'^ 

a + 3H,+ -^--p. (97) 

The effective scalar potential is plotted in Fig. [lO] for A = 1, and it is multivalued indeed. It has a minimum 
depending on the values of n and A. For generic values of the parameters, the minimum of the potential is close to 
(7 = 0, corresponding to infinitely large curvature R = ±oo. Thus, while the field is evolving towards the minimum, it 
evolves oscillating towards a singular point. We have a stable de Sitter minimum and an unstable de Sitter maximum. 
The point R ^ corresponds to a fiat spacetime, which, although a solution for the model, is unstable. Cusps that 
occur at i? = ±l/\/3 are critical points with /" — 0. However, depending upon the values of the parameters, we can 
choose a finite range of initial conditions for which scalar field a evolves to the minimum of the potential without 
hitting the singularity. Note again, that as it was demonstrated in second section, it is just manifestation of II Type 
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finite-time singularity. In Fig. [TTl the cosmological behavior of a is plotted. It is clear the contribution of matter in 
changing the local behavior. 

The time-time component of the equation of motion ^ gives the Hubble equation 

, d{\nF'{R)) I F{R)-F{R)R _ k\ 
^ + It ^ + 6 J{R) (^^^ 

The Einstein gravity is recovered in the limit /' = 1 . The picture of dynamics which appears here is the following: 
above infrared scale (i?c), the expansion rate is set by the matter density and once the local curvature falls below Rc 
the expansion rate gets the effects of modified gravity. 

For pressureless dust, the effective potential presents an extremum at 

2F{R)-RF\R)^~'^p. (99) 

For viable late time cosmologies, the field evolves near the minimum of the effective potential. The finite-time 
singularity, which occurs for the class of models under consideration, severely constrains the field dynamics. 



C. Adding higher curvature corrections to cure singularities 



It is well-known that for large curvature regimes, the quantum effects become important leading to higher curvature 
corrections. Then the program of this paper can be enhanced by considering if higher curvature corrections (as is 
already stressed in third section) added to the original models can solve the singularity problem. In general, higher 
curvature corrections change the structure of the effective potential around the singularity [131 . Keeping this in mind, 
let us consider the modification of the model l(82l) . Although we focus on a specific model, similar results hold for the 
models of the classes considered here. In cosmology, higher curvature corrections appear as Lagrangian contributions 
like L — aiR? -I- a^^R? 4- • • • , and so the most natural choice for the leading order term is aR? /R^} It is well-known 
that the R^ - term may be responsible for infiation in the early Universe if Rc is set to be at infiationary scale [6Q| . 
In the case which we are considering, we have 

nR) = - ',^, +^^, + ^. (100) 



then the field a becomes 



(^)2n + 1 Rl Rc 



, , R 2nX 2nX 
a(R)^2 — a 5 + ^ 101 

^R-J £(^(^)2n + lj 



Higher order corrections naturally include terms like RfivR'^" , but in this paper we focus on the /(i?)-type modified gravity and hence 
simply assume that the corrections are also given by a function of the Ricci scalar. 
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Figure 12: Parametric plot of the effective potential for n = 2, A = 2 and a = 1/2 in presence of correction in two and 
three dimensions for the model l[82|l . 

When \R\ is large in modulus, the first term which comes from ai?^ dominates. In this case, the curvature singularity 
R = ±oo corresponds to cr = ±oo. Hence, by this modification, the minimum of the effective potential is separated 
from the curvature singularity by the infinite distance in the {cr, V{a)} plane. 

For n = 2, X = 2 and a — 0.5, we have a large range of the initial condition for which the scalar field evolves to the 
minimum of the potential as shown in Fig. [I2l In conclusion, the introduction of term formally allows to avoid 
the singularity as it was suggested earlier also in [stI. lisl. fi^ . 

V. CONCLUSIONS 

In this paper we have discussed the future finite-time singularities which can, in principle, appear in dark energy 
universe coming from modified gravity as well as in other dark energy theories. Considering f{R) gravity models that 
satisfy cosmological viability conditions (chameleon mechanism), it is possible to show that finite-time singularities 
emerge in several cases. Such singularities can be classified according to the values of the scale factor a(i), the 
density p and the pressure p as done in SectUll To avoid the singularities, suitable boundary conditions have to be 
imposed which depend, in general, from the parameters of the model as seen in Sect lllll It is interesting that in 
static spherically-symmetric spacetime the finite-time singularity manifests itself as singularity at some specific value 
of curvature. 

Besides, the problem can be analyzed by considering the mass of an auxiliary scalar field coming from the further 
degrees of freedom of /(i?) -gravity. Such a scalar field is heavy in the high-curvature regime whose density is much 
larger than the present cosmological density. Such a field allows to study the singularity problem using the conformal 
transformations. In this case, we have to consider singularities of the scalar field and the related effective potential 
and try to see if they can be avoided in the conformal picture. For example, the most striking feature of the potential 
in Figure m and the core of the problem for infrared-modified f{R) models, is that curvature singularity is at a finite 
distance both in field and energy values. The scalar field a directly feels the matter distribution; for suitable values 
of the parameters, the force is directed to the right, and drives the field a up the wall toward infinite curvature. The 
characteristic scale of the potential V{a) is the curvature Rc which is of the same order of magnitude of the today 
observed cosmological constant. Such a value is extremely low compared to the matter densities we encounter at local 
Solar System and Galactic scales. Given the scales involved, it is easy to drive the scalar field in order to jump the 
potential well considering the dynamics of standard matter which could cause catastrophic curvature singularity. If 
this were to happen, this is in contrast with any viable model. Similarly, matter with sufficiently stiff equation of state 
can destabilize the model by driving the field to unstable points. It is also remarkable that even in fiat spacetime, 
some classes of the classical potential may bring the theory to finite-time singularity. 

Finally we have discussed the possibility to address and solve the singularity problem adopting higher-curvature 
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corrections. In such a way, the future cosmological era has no singularity as it was demonstrated earUer in [stI. [isl. li^l . 
Then, it also does not manifest itself in spherically-symmetric spacetime. The interesting result of this analysis is the 
fact that finite-time singularities, which present in some modified gravity, may not influence the conjectured problem 
with relativistic star formation indicated in ref.fi^l. Of course, this fact depends on the specific characteristics of the 
adopted modified gravity model which can be even totally free of future, finite-time singularities. 
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